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ABSTRACT 

This  paper  compares  approaches  proposed  in  the  literature  for  the 

Capacitated  Plant  Location  Problem.  The  comparison  is  based  on  new 

theoretical  and  computational  results.  The  main  emphasis  of  the  paper  is  on 

relaxations.  In  particular  we  identify  dominance  relations  among  the  various 

relaxations  found  in  the  literature.  We  also  perform  a  probabilistic  study  of 

these  relaxations  using  a  Euclidean  model.  In  the  computational  study,  we 

compare  the  relaxations  as  a  function  of  various  characteristics  of  the  test 

problems.  Several  of  these  relaxations  can  be  used  to  generate  heuristic 

feasible  solutions  that  are  better  than  the  classical  greedy  or  interchange- 

type  heuristics,  both  in  computing  time  and  in  the  quality  of  the  solutions 
found.  /  r-tejer  J,  r.^r  jpr°^ 
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1.  Introduction 

For  many  organizations  the  location  of  plants,  such  as  warehouses  or 
factories,  affects  the  cost  of  supplying  commodities  to  clients  through  the 
transportation  costs  on  the  one  hand  and  the  fixed  costs  of  opening  and 
operating  the  plants  on  the  other  hand.  These  costs  vary  with  the  location 
and  size  of  the  plants.  In  this  paper  we  study  a  classical  model  in  location 
theory. 

The  Capacitated  Plant  Location  Problem  (CPLP)  can  be  formulated  as  the 
mixed  integer  linear  program 


Z  =  Min  Z.Z  .  c. .x. .  +  Z  .  f .y . 

i  J  ij  1J  J  J  J 


Z  .  x.  .  =  1  for  all  i 

J  ij 


(C)  z.  d.x..  <  s  .y  .  for  all  j 

i  i  ij  J  J 


(N)  0  <  x _  <  1 ,  0  <  y^  <  1  for  all  i  and  j 


(I)  y  integer  for  all  j, 


where  d^^  is  the  demand  of  client  i,  fj  and  Sj  are  the  operating  cost  and 
capacity  of  facility  j,  if  it  is  open,  and  c^j/d^  is  the  unit  transportation 
cost  between  i  and  j.  The  decisions  to  be  made  are  represented  by  the 
variables  yj  and  x^j,  namely  y j  =  1  if  facility  j  is  open  and  0  otherwise,  and 
x^  is  the  fraction  of  the  demand  of  client  i  met  from  facility  j.  The 
constraints  (D)  are  the  demand  constraints,  (C)  the  capacity  constraints,  (N) 


the  nonnegativity  and  simple  upper  bound  constraints  and  (I)  the  integrality 
constraints . 

The  applications  of  this  model  are  not  limited  to  plant  location.  For 
example,  the  same  mathematical  model  is  appropriate  for  making  optimal  lot 
sizing  decisions  in  production  planning.  If  d^  is  the  demand  in  period  i,  f j 
and  s^  are  the  set  up  cost  and  production  capacity  in  period  j,  then  CPLP  is  a 
valid  model  where  the  decision  variables  are  interpreted  as  follows:  x^j  is 
the  fraction  of  the  demand  d^  produced  in  period  j  and  yj  =  1  if  production 
occurs  in  period  j,  0  otherwise.  Finally,  in  this  model,  c^j/d^  is  the  unit 
production  cost  in  period  j  plus  the  inventory  holding  cost  until  period  i,  if 
i  >  j  (or  the  backlogging  cost  if  i  <  j).  The  model  CPLP  has  also  been 
proposed  in  the  context  of  telecommunication  network  design  (Kochman  and 
McCallum  [22]  and  Mirzaian  [24]),  machine  replacement  and  vehicle  routing  (see 
Sridharan  [28]). 

The  formulation  of  the  CPLP  can  be  enriched  by  simple  valid 
inequalities.  For  example,  the  variable  upper  bound  constraints 

(B)  x.  .  <  y .  for  all  i, 1 

1 J  J 

and  the  total  demand  constraint 

(T)  Ej  SjYj  2  <V 

To  see  that  (B)  is  valid  for  CPLP,  note  that  x^j  <  1  by  constraint  (N).  Using 

(C)  and  (I)  it  follows  that  x^j  <  y j .  To  see  that  (T)  is  valid,  sum  up  the 
constraints  (C)  and  use  (D). 


One  obtains  various  lower  bounds  on  Z  by  relaxing  subsets  of  the 
constraints  (D),  (C),  (I),  (B)  and  (T)  either  completely  or  in  a  Lagrangian 
fashion.  We  use  the  following  notation.  If  a  set  of  constraints,  say  (A),  is 
relaxed  completely  the  resulting  bound  is  denoted  by  Z^.  If  constraints  other 
than  (I)  are  relaxed  in  a  Lagrangian  fashion,  the  corresponding  Lagrangian 
dual  bound  is  denoted  by  Z^.  For  example 


Zn  =  Max  Min  Z .  Z  .  c..x..  +  Z  .  f.y.  -  Z .  u.(Z.  x..  -1) 
D  u  x,y  1  J  LJ  iJ  J  J  J  1  1  J  iJ 


Z .  d . x .  .  <  s  y  .  for  all  j 
i  i  ij  J  J 


0  <  x .  .  <  y  .  for  all  i , j 

zj 


Z  .  s  ,y  .  >  Z.d. 
j  rj  z  z 


y^  e  {0,1}  for  all  j . 


These  lower  bounds  provide  the  basis  for  numerous  algorithms.  For  example, 
the  following  relaxations  have  been  used  in  branch  and  bound  algorithms. 

Z  (This  relaxation  is  the  so-called  weak  linear  programming  relaxation 
of  CPLP.  It  was  used  by  Sa  [27]  and  improved  by  Akinc  and  Khumawala 
[2]  using  ad  hoc  rules). 

(This  is  the  strong  linear  programming  relaxation  of  the  problem.  It 
was  used  by  Davis  and  Ray  [7],  Guignard  and  Spielberg  [18],  and  Baker 


ZTp  (Geoffrion  and  Me  Bride  [14]) 
Zq  (Nauss[25]) 
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Z'q  (Christof ides  and  Beasley  [6].  In  this  relaxation,  the  total  demand 
constraint  (T)  is  replaced  by  z  y  >  K  where  K  is  the  smallest 
integer  such  that  the  sum  of  the  K  largest  Sj's  is  at  least  Z^d^). 

Zq  (Van  Roy  [301  :  Guignard  and  Kim  [16]). 

Recently,  a  new  type  of  relaxation  has  been  introduced  by  several 
authors.  The  method,  known  as  variable  splitting,  can  be  defined  in  general 
for  (mixed)  integer  programs.  It  is  obtained  by  rewriting  the  program 

Z  =  min  c  u 

(51 )  A1  u  <  b1 

( 52 )  A2  u  ^  b2 

u  e  X 
as 

Z  =  min  c  u 

A^  u  <  b^ 
u  e  X 
v  -  t>2 
u  e  X 
u  =  u 

and  by  performing  a  Lagrangian  relaxation  of  the  equalities  u  =  u.  The 
Lagrangian  dual  bound  obtained  in  this  way  has  the  value 

Z  (S^ ,  S2)  =  min  c  u 

ue  (A^u  <  b^,  ueX)  n  (A^u  <  b^ ,  ueX) 


whereas  the  classical  Lagrangian  bound  obtained  by  dualizing  (S^)  is 


Z  =  min  cu 
b1 

ue  {A^u  <  b^}  n  (A^u  <  b^,  ueX} 

(see  Guignard  and  Kim  [17]  and  Geoffrion  [13]  respectively).  As  a 
consequence,  <  Z(S^,S._,),  i.e.  the  bound  obtained  by  variable  splitting 
is  always  as  tight  as  the  corresponding  Lagrangian  bound. 

In  a  recent  paper,  Barcelo,  Fernandez  and  Jonsten  [4]  propose  the 
following  split  for  CPLP.  (S^)  corresponds  to  the  demand  constraints  (D),  X 
to  the  nonnegativity  and  0,  1  requirements  (N)  and  (I),  and  (S2)  to  the  other 
constraints,  namely  (C),  (B)  and  (T).  We  show  in  Section  2  that  the  proposed 
relaxation  yields  the  bound  Zj>.  Other  splits  yield  other  relaxations. 

Our  motivation  for  this  study  is  twofold,  theoretical  and  computa¬ 
tional.  On  the  theoretical  side,  we  rank  relaxations  according  to  the 
tightness  of  the  bounds  that  they  yield.  This  analysis  is  performed  in 
section  2.  Several  inequalities  given  in  that  section  are  classical,  but  we 
have  had  a  few  surprises.  For  example,  consider  Zq,  Z^q  and  Z^-,.  It  is  easy 
to  show  that  >  Z^  >  zj  and  both  Van  Roy  [30]  and  Guignard  and  Kim  [16] 
have  noticed  that  the  bound  Zq  is  often  strictly  stronger  that  Z^  in 
practice.  However  Zq  is  also  harder  to  compute  than  Z ^  and,  in  both 
references,  the  Lagrangian  relaxation  of  (T)  is  introduced.  Indeed,  it  seems 
that  computing  Z^q  is  a  good  compromise  between  bound  quality  and  ease  of 
solution.  Our  Theorem  of  Section  2  shows  that,  in  fact,  Z^q  =  Zq.  So  Z^q 


til* 


cannoc  be  used  as  a  shortcut  in  approaching  the  quality  of  the  bound  Zq.  On 
the  other  hand,  the  solution  of  by  subgradient  optimization  may  still  have 
merits  as  it  is  possible  that  it  converges  faster  to  the  common  bound  =  Zq 
than  a  subgradient  approach  applied  directly  to  Z£. 

Another  interesting  result  from  Section  2  is  that  variable  splitting  does 
not  yield  stronger  bounds  than  Zq. 

Section  3  analyzes  the  complexity  of  solving  the  different  Lagrangian 
relaxations.  We  also  consider  the  complexity  of  the  problem  CPLP  itself  under 
two.  different  assumptions  on  the  capacities  of  the  facilities. 

Section  4  contains  our  computational  experience  with  various 
relaxations.  Nearly  all  the  researchers  in  the  field  have  tested  their 
algorithms  on  the  same  set  of  data,  the  famous  Kuehn  and  Hamburger  test 
problems  as  adapted  by  Khumawala  [21].  These  test  problems  have  been  a  great 
help  in  comparing  the  various  solution  methods.  At  the  same  time,  it  should 
be  noted  that  this  test  set  is  relatively  small  and  has  a  very  special 
property  :  all  the  facilities  have  the  same  capacity  and  operating  cost. 
Therefore  it  seems  appropriate  to  conduct  a  computational  study  on  a  different 
and  larger  set  of  test  problems.  These  problems  have  varying  capacities  and 
operating  costs  but,  for  each  facility,  the  operating  cost  is  correlated  to 
the  capacity  according  to  an  increasing  concave  function  (economies  of 
scale).  The  duality  gap  yielded  by  various  relaxations  appears  to  be 
sensitive  to  characteristics  of  the  data  such  as  total  capacity  versus  total 
demand.  To  demonstrate  such  relationships  we  perform  our  computational 


experiments  as  a  function  of  o  =  I^s^/E^d^.  Keeping  everything  else  equal, 
we  scale  the  Sj's  in  order  to  get  different  values  of  o.  For  example  we  find 
that  adding  the  total  demand  constraint  (T)  improves  the  relaxations 
significantly  when  p  is  small. 


Our  computational  experience  also  covers  some  existing  heuristics.  Here 
as  well  we  have  had  surprises.  A  number  of  heuristics  which  perform  well  on 
Kuehn  and  Hamburger  problems  have  a  poor  performance  on  our  test  problems.  We 
recommend  the  use  of  relaxation-based  heuristics  rather  than  the  usual  greedy 
or  interchange-type  heuristics.  Consider  any  iterative  method  (such  as 
subgradient  optimization)  to  compute  Z^.  For  each  ufc,  t  >  1,  generated  during 
the  iterative  process,  the  subproblem  in  x  and  y 


Min  Z.Z.  c..x..  +  Z.  f  .y  .  -  Z.  u.(Z.  x..  -1) 
x,y  i  J  ij  ij  J  J  J  i  i  J  i-J 


Z .  d .  x .  .  <  s  .y  .  for  all  j 

i  i  ij  r  j 


0  <  x.  .  <  y  . 
ij 


for  all  i,  j 


Z  .  s  .y  .  >  Z.d. 
J  J  J  i-i 


y.  t  {0,1} 


for  all  j 


is  solved.  The  vector  (x^,  y q )  so  obtained  produces  a  Lagrangian  bound,  but 

the  vector  yD  can  also  be  used  to  generate  a  heuristic  solution.  Indeed, 

given  y^,  one  can  solve  a  transportation  problem  to  get  a  feasible  solution 

Xp  .  Thus  at  each  iteration  of  the  iterative  process,  a  heuristic  solution 

(xp,  y^)  can  be  generated.  The  best  solution  found  in  50  or  so  subgradient 

iterations  can  be  used  as  a  heuristic.  We  find  that  this  heuristic  is  just  as 

fast  and  gives  much  better  solutions  than  the  more  commonly  used  greedy- 

T 

interchange  approach.  Instead  of  Zp,  one  can  also  use  Z^  or  Z^  as  the  basis 
for  the  heuristic  but  we  do  not  recommend  it,  as  little  additional  accuracy  is 


gained  at  much  computational  cost. 


Section  5  contains  probabilistic  results.  This  analysis  improves  our 
understanding  of  the  CPLP  in  a  new  dimension  and  complements  the  computational 
study  of  Section  4.  For  the  probabilistic  mode1  that  we  choose,  the 
relaxations  fall  into  three  categories.  For  relaxations  such  as'  Z  ,  the 
relative  error  as  an  estimate  of  Z  can  be  expected  to  be  over  66%.  For  bounds 
such  as  Z *  or  Zp,  the  relative  error  can  be  expected  to  be  below  0.2%. 
Finally  for  the  relaxations  Zq,  Z ^  and  Z£,  the  relative  error  can  be  expected 
to  go  to  0  as  the  size  of  the  instances  goes  to  infinity.  In  particular,  we 
show  that  the  classical  linear  programming  relaxation  Z1  has  a  "duality  gap" 
which  is  no  greater  for  CPLP  than  for  the  Simple  Plant  Location  Problem 
(SPLP).  This  seems  to  contradict  the  common  belief  that  CPLP's  have  larger 
duality  gaps.  We  should  remember  that  our  probabilistic  analysis  is  limited 
by  various  assumptions.  For  example,  it  assumes  that  the  number  of  clients 
served  by  each  open  facility  is  large. 

2.  Relative  Strength  of  the  Relaxations 

In  Theorem  1,  we  compare  the  strength  of  the  various  linear  programming 
relaxations  and  lagrangian  duals  of  CPLP.  We  compare  all  the  possible  bounds 
that  can  arise  in  this  way,  except  that  we  do  not  relax  (N)  (there  is  never 
any  computational  gain  in  doing  so),  nor  do  we  consider  complete  relaxations 
of  (C)  or  (D)  since  we  feel  that  such  bounds  would  be  of  little  interest  in 
solving  CPLP  (ZC  and  Z®  are  NP-hard  to  compute,  yet  they  lack  important 
aspects  of  CPLP.)  Specifically,  we  consider  all  the  combinations  where  (T) 
and  (B)  are  present,  relaxed  completely  or  in  a  Lagrangian  fashion,  and  where 
(C)  and  (D)  are  present  or  relaxed  in  a  Lagrangian  fashion,  i.e.  3x3 x2x2  =  36 
relaxations,  as  well  as  the  LP  relaxations  Z*,  Z®^,  Z^  and  z"®®* .  In 


addition,  we  consider  the  relaxation  Zp  of  Christofides  and  Beasley  defined  in 


ir 


Section  1  above.  Theorem  1  states  that  these  41  relaxations  only  yield  7 
genuinely  different  bounds.  It  is  interesting  to  note  that  6  of  these  7 
bounds  have  already  all  been  proposed  in  the  literature  as  the  basis  for 
branch-and-bound  algorithms. 

In  Theorem  2,  we  briefly  discuss  relaxations  that  can  be  obtained  by 
variable  splitting.  We  show  that  no  new  bound  arises  from  these  relaxations. 

We  use  the  following  notation.  Given  a  constraint  set  (S),  P(S)  denotes 
the  convex  hull  of  the  vectors  satisfying  (S).  More  generally,  if 
(S^ ),...,  (S^)  are  constraint  sets,  P(S1...S,K)  denotes  the  convex  hull  of  the 
vectors  satisfying  (S1 ),..., (S^) .  For  example,  with  this  notation,  we  can 
write 


~TB  .  „ 

Z  =  mm  cx  +  fy 

( x ,y )  e  P(DCIN). 

Using  a  classical  theorem  (see  Geoffrion  [13,  Theorem  ID]  for  example), 
the  Lagrangian  duals  can  also  be  expressed  using  this  notation.  As  an 
example, 

ZTC  =  min  cx  +  fy 

(x,y)  £  P(TC )  n  P(DINB) . 

Now  we  are  ready  to  prove  the  following  Theorem. 

Theorem  1  The  following  inequalities  hold 


y.viY- 


-V-  »  -Vl.--  - 
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(2.1) 


,BI 


SZ^  Z'  <1zdS]zc<z, 


every  inequality  in  (2.1)  is  strict  for  at  least  one  instance  of  CPLP,  and  no 
other  inequality  exists  among  the  values  of  (2.1)  except  those  derived  by 
transitivity . 

The  following  equalities  hold 


(2.2)  1  ,  ZB  =  ZB  =  ZT  =  ZT  =  ZTB  =  Z®  =  ZTb  ZTB 


(2.3) 


7  -  7  -7  -7  -  Z 

D  DC  '  BD  "  BC  8DC 


=  Z 


(2.4) 


ZB  _  2® 

C  ‘  DC 


(2.5) 


zi  =  Z 


TC 


(2.6)  Z1  =  ZTI 


T  T  T 

ZTD  =  ZD  =  ZTDC  =  ZDC  =  ZBTC  :  ZBC  "  ^BTDC 


'BDC 


7b  7bt  7  7t 

'  TD  ■  D  "  BTD  "  BD 


(2.7) 


BI  TBI  B  yTB  B  TB 

TC  C  *  TDC  ‘  DC' 


Proof :  In  (2.1),  the  inequalities  Z®1  <  Z^  and  z£  <  2^  <  Z  are  immediate  by 
relaxation.  Z^  <  Z^  follows  from  the  fact  that  the  constraint  used  in  Z'p 
in  place  of  (T)  is  actually  implied  by  (I)  and  (T).  Using  the  property  that 
Lagrangian  duals  are  at  least  as  strong  as  the  corresponding  linear 


and 


programming  relaxation  (Geoffrion  [13]),  we  get  Z  <  Z^,  Z  <  Z^ 


IT  T 

Z  <  Zr. 


The 


last  inequality  implies  Z1  <  Z^  provided  we  can 


show 


Z^  =  ZIT.  This  will  be  done  in  (2.6) 


,B 


Z^  <  Zp  follows  from 


B  B  B 

<  Zn  provided  we  can  show  lr  -  Znr.  This  will  be  done  in  (2.4). 

Z_.  This 


'DC  "  "D  - -  ""  -  “C  "  "DC 

Finally  Z <  Z c  follows  from  ZDC  <  Z c  provided  we  can  show  ZD  -  ~DC . 

will  be  done  in  (2.3).  So,  to  show  the  validity  of  the  theorem,  it  remains  to 

prove  the  equalities  (2.2)  -  (2.7)  and  to  exhibit  instances  showing  that  all 

the  inequalities  in  (2.1)  can  be  strict  and  that  no  other  inequality  exists 

among  the  bounds  considered  in  (2.1).  These  results  are  provided  in  the  next 

two  lemmas .  □ 


Lemma  1  The  equalities  (2.2)  -  (2.7)  hold. 

Proof :  We  give  polyhedral  proofs  of  these  results.  See  [28]  for  algebraic 
proofs . 

The  constraints  (T)  and  (B)  can  be  derived  from  (D),  (C),  (I)  and  (N)  as 
shown  in  the  introduction.  Therefore  P(DCIN)  =  P(TBDCIN),  proving  (2.2). 

To  prove  (2.3),  we  first  show  that  =  ^BDC  ancl  tlien  ZD  =  ZD‘ 

The  first  equality  holds  if  we  can  prove  the  polyhedral  result 

(2.8)  P(D)  n  P(BCINT)  =  P(BDC)  n  P(INT). 

Since  P(BDC)  =  P(D)  n  P(BC)  we  only  have  to  show 

Claim 


(2.9) 


P(BC)  n  P(INT)  =  P(BCINT) . 


Proof  of  Claim: 


By  definition  of  P,  P(BCINT)  C  P(BC)  n  P(INT).  Now  consider  any  extreme  point 

of  P(BC)  n  P(INT),  say  (x,y).  If  y  has  only  0,1  coordinates,  then 

(x,y)eP(BCINT)  and  the  claim  is  proved.  So  assume  that  y  has  at  least  one 

fractional  component.  We  will  show  that  this  cannot  occur.  Denote  by 

Py(INT)  the  projection  of  the  polyhedron  P(INT)  onto  the  y-space.  Since 

Py(INT)  has  only  integral  vertices,  y  is  not  a  vertex  and  therefore 

y  =^y1  +  |y2  where  y',y2  e  P  (INT).  If  y.  -  0,  set  x^  -  x^\  =  0. 

Else,  set  x^ .  =  x.  .  (y^/y.)  for  k=1,2.  We  have  (x^,yk)  e  P(BC)  n  P(INT) 

ij  ij  J  J 

-  —  1  11  12  2 

for  k= 1 , 2  and  (x,y)  =  -  (x  ,y  )  +  -(x  ,y  ).  This  contradicts  the  assumption 
that  (x,y)  is  an  extreme  point  of  P(BC)  n  P(INT).  Thus  the  claim  is 
proved. 

g 

Now  we  prove  ZQ  =  Z D  by  showing 

(2.10)  P(D)  n  P(BCINT)  =  P(D)  n  P(CINT). 

This  relation  follows  from  the  fact  that  (C),  (I)  and  (N)  imply  (B).  This 
completes  the  proof  of  (2.3). 

To  prove  (2.4),  we  show  that 

(2.11)  P(C)  n  P(INTD)  =  P(DC)  n  P(INT). 

Since  P(DC)  =  P(C)  n  P(D),  it  suffices  to  show  P(D)  n  P(INT)  =  P(DINT).  The 
constraints  (D),  (I),  (N)  and  (T)  completely  separate  into  the  constraints 
involving  x  and  those  involving  y.  The  convexification  step  only  involves  the 


variables  y  and,  therefore,  the  result  follows.  This  completes  the  proof  of 

(2.4). 

The  equality  (2.5)  follows  from  P(C)  n  P(DINB)  =  P(TC)  n  P(DINB),  since 
(T)  is  a  consequence  of  (C)  and  (D). 

Mow  we  show  (2.6).  First,  Z^  =  Z since  (T)  follows  from  (C)  and  (D). 

BT 

First,  we  give  a  polyhedral  proof  of  Z^  =  Z^^.  We  have 

P(D)  n  P(CIN)  =  P(D)  n  P(CINB)  =  P(TD)  n  P(CIMB),  where  the  first  equality 
holds  because  (C),  (I),  (M)  imply  (B)  and  the  second  equality  holds  because 

(D)  and  (C)  imply  (T).  So  all  the  relaxations  in  (2.6)  are  at  least  as  strong 

ITT  T 

as  Z  =  Z  To  complete  the  proof  of  (2.6)  if  suffices  to  show  Z^^  =  Z 

and  ~  Z  * 

To  show  Z^q  =  Z1,  note  that 

(2.12)  P(TD)  n  P(CINB)  =  P(TDCMB) 

follows  from  P(CINB)  =  P(CMB).  To  show  Zg^Q  -  Z*,  note  that 

(2.13)  P(BTC)  n  P(DIN)  =  P(BTCDN) 

follows  from  P(DIN)  =  P(DN). 

Finally  consider  (2.7).  Z®^  =  Z^®*  holds  as  (T)  follows  from  (C)  and 
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(D).  So  it  suffices  to  show  Z  =  2°^.  Similar  to  the  proof  of  (2.13), 

this  follows  from  P(DIN)  =  P(DN).  Q 

To  complete  the  proof  of  the  theorem,  it  suffices  to  exhibit  instances 


showing  that  all  the  inequalities  can  be  strict. 


Lemma  2 :  There  exist  instances  of  the  CPLP  showing  that  each  of  the  following 
inequalities  can  hold. 


(2.14)  ZBI  <  Z1 

(2.15)  <  ZB 

(2.16)  zj  <  ZB 

(2.17)  zj  <  Zp 

(2.18)  ZQ  <  Z7C 

(2.19)  ZB  <  Z1 

(2.20)  Zc  <  Z. 


Proof :  First  we  give  an  example  that  satisfies  (2.14)  and  (2.15). 

1 

s  =  3  5 
f  =  2  3 

BI 

To  compute  Z  ,  we  solve  the  corresponding  linear  program.  Its  optimum 
solution  is  y1  =  1 ,  y2  =  2/5,  xn  =  1,  x12  -  0,  x21  =  1/3  and  x22  -  2/3, 
yielding  the  bound  Z®*  =  16/5.  Now  consider  Z*.  Its  optimum  solution  is  the 
same  as  for  Z®^ ,  except  that  y2  =  2/3.  This  yields  Z^  =  4.  Thus  this  example 
shows  that  (2.14)  can  occur. 

Note  that,  for  this  instance,  the  constraint  (T')  used  in  Zq  is  y1  +  y2  > 
1  and  therefore  does  not  strengthen  the  relaxation  Zp.  Since  we  have  proved 
that  Zq  =  Z1,  we  get  Z'D=Z2=ZI=4. 

g 

Now  we  compute  Z^,.  Consider  the  Lagrangian  relaxation 


=  Min  z.l.  c.  .x.  .  +  L  .  f.y.  +  z.v.U.d.x..  -  s.y.' 
Cv  iJ  iJiJ  J  J  J  J  J  i  i  ij  /j 
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l.x.  . 
j  ij 


=  1 


z . 
J 


s  .y  .  >  z . 
J  J  1 

e  {0,1} 


d. 

i 


for  all  i 


for  all  j . 


If  we  set  v1  =  2/3,  v2  =  0,  an  optimum  solution  of  the  Lagrangian 


relaxation  is  =  y2  =  1,  x^  =  1,  x12  =  0,  x21  =  1/3,  x22  =  2/3.  The 


corresponding  value  is  =  5.  Since  we  also  have  Z  =  5,  we  deduce 


Z^  =  Z  =  5.  Therefore  the  strict  inequality  (2.15)  is  satisfied  by  this 


instance  of  the  CPLP. 

To  prove  that  (2.16)  and  (2.17)  can  hold,  consider  the  previous  example, 
where  s2  =  5  is  changed  into  s2  =4.  It  is  easy  to  check  that  this 


D  r 

modification  does  not  affect  Z^  and  Z.  But,  now,  the  constraint  ( T  * )  used 


B 


in  Zp  is  y1  +  y2  >  2  and  therefore  Z^  =  Z^  =  5.  Next  we  compute 


zc  5  Ma\>o  zcv’  “here 


zj,  =  Min  Z.Z.  c.  .x.  .  +  I.f.y.  +  I.v.(z.d.x..  -  s.y.) 
Cv  l  j  ij  ij  J  JJJ  j  j  l  l  ij  JJJ 


z  .  x .  .  =  1  for  all  i 
J  ij 

0  <  x .  .  <  y .  for  all  i ,  j 
ij  J 

y^  e  {0,1}  for  all  j. 

This  constraint  set  is  totally  unimodular  as  it  is  a  simple  plant  location 

problem  with  two  clients  and  two  locations  (see  Cho,  Padberg,  Johnson  and 

T  I 

Rao[5]  for  example).  Therefore  Z^  =  Z  =4.  This  proves  that  (2.16)  and 
(2.17)  can  occur. 

Next  we  give  an  example  that  satisfies  (2.18). 

d- 

10  0  1 

(c  )  =010  1 

0  0  1  1 
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s  =  3  3  3 

f  =  1  1  1 

The  optimum  value  of  the  linear  programming  relaxation  Z1  -  3/2  is  obtained  by 

taking  the  solution  y  j  =  1/2,  Xjj  =  1/2  if  i  *  j  and  x.j  =  0  if  i  =  j. 

Note  that  for  this  instance  Z^j  =  Zq  =  Zp  since  there  is  no  capacity 

T  I 

constraint  in  effect.  Since  we  have  proved  that  Z  q  =  Z  ,  we  obtain  Z^  =  Z^  = 
ZTd  =  Z1  =  3/2. 

Similarly  Z^  -  Z  since  the  capacity  constraints  are  automatically 

satisfied.  This  shows  that  Z^  =  Zq  =  z£  =  Z  =  2. 

Therefore  this  instance  shows  that  possibility  (2.18)  can  occur. 

Finally  we  give  an  instance  that  satisfies  (2.19)  and  (2.20). 

d= 

0  2  1  2 

(c.j)  =102  2 

2  10  2 


s  =  3  3  3 

f  =  3/2  3/2  3/2 


8  2 

To  prove  (2.19),  note  that  Z^  =  3  since  y,=  y^=  y^  =  j  anci  xij  =  1  ^  =  J  anci  0 

if  i*j  belongs  to  P(CD)  n  P(INT).  On  the  other  hand,  Z*=4  since 

y1"  y2=  y3=  3  and  x11=  X22=  X33=  3’  X13=  X21=  x32=  3’  X12  =  X23~  X31~  °  1S  an 
optimum  solution  of  the  linear  programming  relaxation.  To  prove  (2.20),  note 

that  Z  =  9/2  is  obtained  by  setting  y1  =  y£  =  y^  =  1,  x^j  =  1  if  i  =  j  and  0 

if  i  *  j.  For  the  solution  of  Zq,  consider  v1  -  ^2  =  v3  '  1/6  in  the 

Lagrangian  relaxation  Zqv.  Then  the  same  solution  (x,y)  is  optimum.  It 

yields  the  value  Z^v  =  4.  Furthermore  the  subgradient  at  this  point  is 

negative  in  each  component,  so  its  projection  on  the  nonnegative  orthant  is 

zero.  This  shows  that  Zq  =  4,  proving  (2.20). 


» 

7. 

£ 

/ 


IO 

This  completes  the  proof  of  Lemma  2.  □ 

We  close  this  section  with  results  about  variable  splitting.  It  will  be 

useful  to  generalize  the  notation  introduced  in  Section  1.  If  the  set  of 

1  k 

inequalities  (S^)  consists  of  the  constraint  sets  (S^ (S^ )  and  (S2) 

1  z 

consists  of  (S^) , . . . , (S^) ,  then  the  bound  defined  by  the  split  ( S ^ ) ,  (S2) 

1  k  1  S. 

and  X  is  denoted  by  2(S^...S^,  S^...S^).  As  mentioned  in  the  introduction, 
this  value  is  obtained  by  minimizing  ex  +  fy  over 

(2.21)  (x,y)  e  P(sj...S^  X)  n  PtS^.-.S*  X). 

Let  us  consider  all  the  possible  splits  of  the  constraints  (D),  (C),  (B),  (T), 
(I)  and  (N)  for  CPLP.  Our  main  result  is  that  no  new  bound  arises  this  way. 
Therefore  variable  splitting  can  be  viewed  as  a  different  method  for  obtaining 
the  bounds  presented  in  Theorem  1.  We  make  the  following  observations.  If 
(I)  is  not  in  X,  then  only  one  of  the  two  polyhedra  in  (2.21)  is  convexified 


by  integality  and  therefore  the  corresponding  bound  is  one  of  the  Lagrangian 
dual  bounds  studied  in  Theorem  1.  So  we  assume  that  (I)  belongs  to  X.  We 
also  assume  that  (N)  belongs  to  X  as  no  benefit  can  be  gained  by  doing 
otherwise.  There  are  25  possible  splits  of  (D),  (C),  (B),  (T),  (I)  and  (N) 
with  the  property  that  (I)  and  (N)  belong  to  X.  If  ( S 1 )  or  (S2)  is  reduced  to 
one  of  the  sets  (B),  (T)  or  (B)  U  (T),  then  we  do  not  get  a  true  relaxation  of 
CPLP,  as  one  of  the  subproblems  is  CPLP  itself.  Therefore  the  true 
relaxations  have  the  property  that  (D)  belongs  to  (S^  and  (C)  belongs  to 
(S2)  or  vice-versa.  There  are  9  such  relaxations  of  CPLP.  The  bounds  they 
produce  equal  Zq  or  Zp. 


Theorem  2  The  following  equalities  hold 


(2.22)  Z(T,»)  =  Z(B,*)  =  Z ( BT , * )  =  Z 

where  *  stands  for  any  subset  of  the  remaining  constraints, 

(2.23)  Z(C,D)  =  Z(C,DB)  =  Z(C,DT)  =  Z(C,DBT)  =  Z(D,TC)  =  Z(BD,TC)  =  Z 

(2.24)  Z ( D , BC )  =  Z(D,TBC)  =  Z(TD,BC)  =  Z 

Proof  =  (2.22)  is  immediate  as  one  of  the  two  subproblems  is  CPLP  itself. 

Next,  we  prove  (2.23).  Note  that  Z(C,D),  Z(C,BD),  Z(C,DT)  and  Z(C,DBT) 
are  all  at  least  as  strong  as  Z c,  and  Z(C,D)  is  the  strongest  of  the  four.  To 
prove  Z^  =  Z(C,D)  we  show  the  polyhedral  result  P(C)  n  P(DINTB)  = 

P(CINTB)  n  P(DINTB) .  This  follows  from  P(CINTB)  =  P(C)  n  P(INTB)  as  we 
have  already  shown  in  (2.9),  and  from  P(DINTB)  C  P(INTB).  To  complete  the 
proof  of  (2.23),  it  remains  to  show  Z(D,TC)  =  Z(BD,TC)  =  Zc>  The  feasible  set 
for  Z(D,TC)  is  P(DINB)  n  P(TCINB).  Using  (2.9),  we  have  P(CINTB)  =  P(B)  n 
P(CINT)  and  therefore  Z(D,TC)  =  Z(BD,TC).  Again,  from  (2.9)  P(DINB)  n 
P(TCINB)  =  P(DINB)  n  P(INT)  n  P(C).  In  order  to  prove  Z(D,TC)  =  Z q,  we  must 
show  P(DINTB)  =  P(DINB)  n  P(INT).  Let  (x,y)  be  an  extreme  point  of  P(DINB) 
n  P(INT).  Note  that  (y:  (x,y)eP(INT)  for  some  x}  C  {y:  (x,y)  e  P(DINB)  for 

some  x}  =  {y:  E j  y j  >  1 }  using  the  fact  that  d^  >  0.  Therefore  y 
an  extreme  point  of  the  knapsack  problem  defined  by  (T),  i.e.  y  is  a  0,1 
vector.  As  (x,y)  e  P(DINB),  x  is  also  a  0,1  vector.  Thus  (x,y)  is  a  0,1 
vector  and  (x,y)  e  P(DINTB). 


!*] 


Finally  consider  (2.24).  Z(D,BC)  and  Z(D,TBC)  are  at  least  as  strong  as 
Zp.  To  show  that  they  are  equal  to  it,  it  suffices  to  show  that  the  feasible 
set  for  Z(D,BC),  namely  P(DINT)  n  P(BCINT),  is  identical  to  the  feasible  set 
for  Zp,  namely  P(D)  n  P(BCINT).  This  follows  from  P(DINT)  =  P(D)  n  P(INT) 

and  P(BCINT)  C  P(INT).  Mow  consider  Z(TD,BC).  It  has  feasible  set 
?(TDIN)  n  P(BCIN).  The  equalities  P(TDIN)  =  P(D)  *  P(INT)  anH 

P(BCIN)  =  P(BC)  n  P( IN)  imply  Z(TC,BC)  =  ZD.  D 

3.  Computational  Complexity 

B  T 

In  Section  2,  seven  different  bounds  were  obtained  for  CPLP,  namely  Z  , 
T  R  T 

Z  ,  Zp,  Zp,  Zq  ,  Zq  and  Zp.  Here  we  consider  the  computational  complexity  of 
calculating  these  bounds.  We  also  consider  the  complexity  of  computing  the 
optimum  value  Z  of  CPLP  itself. 

Theorem  3 

(3.1)  ZBI,  Z1  and  Zp  can  be  computed  in  polynomial  time. 

(3.2)  Zp  and  ZB  can  De  computed  in  pseudo-polynomial  time. 

(3.3)  Computing  Zq,  Zp  or  Z  is  strongly  MP-hard. 

Proof  =  Linear  programs  can  be  solved  in  polynomial  time,  so  (3-1)  holds  for 
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Z  and  Z  .  Zp  is  obtained  by  minimizing  cx  +  fy  over 

P ( D )  n  P(CBINT')  =  P(DCB)  n  P(INT').  The  constraint  set  (N),  (T')  is 
totally  unimodular.  So  computing  Zp  reduces  to  the  solution  of  a  linear 
program.  This  proves  (3.1). 

Computing  Zp  or  Z“  is  NP-hard  since,  when  c=0,  the  problem  is  nothing 


but  a  knapsack  problem.  Next  we  show  that  Zp  can  be  computed  in  pseudo¬ 
polynomial  time.  To  get  Zn,  we  optimize  cx  +  fy  over  P(DCB)  n  P(INT).  If 


follows  from  a  result  of  Grtischel,  Lovasz  and  Schrijver  [15]  that  optimizing 
over  P(DCB)  n  P(INT)  requires  at  most  pseudo-polynomial  time  since 
optimizing  over  P(DCB)  is  polynomial  and  optimizing  over  P(INT)  is  pseudo- 

D 

polynomial.  To  get  1  ,  we  optimize  cx  +  fy  over  P(DC)  n  P(INT).  so  the 
same  argument  can  be  used.  This  proves  (3.2). 

Finally,  (3.3)  follows  from  the  fact  that  SPLP  is  strongly  NP-hard  and 
polynomially  reduces  to  CPLP  by  setting  the  capacities  Sj  large  enough.  □ 

Although  CPLP  is  NP-hard,  there  are  some  special  cases  that  can  be  solved 
in  polynomial  time.  In  the  next  theorem,  we  show  that  CPLP  can  be  solved  in 
polynomial  time  when  the  capacities  and  the  demands  are  taken  from  the  set 
{1,2}  and  it  is  NP-hard  when  the  capacities  are  taken  from  the  set  { 1 , . . . , p } , 
p>  3  even  if  all  the  demands  are  equal  to  1 .  We  also  consider  the  complexity 
of  computing  ZD,  Z£,  Z£  and  Zc  under  the  last  assumption. 

Theorem  4 

(3-4)  CPLP  is  strongly  NP-hard  when  the  capacities  are  taken  from  the  set 

[1,...,p]  for  any  fixed  p  >  3,  even  if  all  the  demand  are  equal  to  1.  Under 

B  T 

this  assumption,  Zq  or  Z£  can  be  computed  in  polynomial  time  but  computing  Zq 

or  Zq  is  strongly  NP-hard. 

(3.5)  CPLP  can  be  solved  in  polynomial  time  when  the  capacities  and 
demands  are  taken  from  the  set  {1,2}. 

Proof :  We  will  prove  that  CPLP  is  strongly  NP-hard  for  p  >  3  by  reducing  the 
3-dimensional  matching  problem  (3DM)  to  CPLP. 

3DM  is  defined  by  three  distinct  sets  I,  J,  K,  each  of  cardinality  n,  and 
a  family  F  of  triplets  (i,j,k)  where  i  e  I,  j  e  J  and  k  e  K.  The  question 


is  whether  there  exists  S  C  F  of  cardinality  n  such  that  the  3n  elements  of  I, 
J,  K  each  appear  exactly  once  in  the  triplets  of  S.  This  problem  is  NP- 
complete  (Garey  and  Johnson  [12]).  We  construct  a  CPLP  as  follows.  There 
are  3n  clients  corresponding  to  the  elements  of  I,  J,  K  and  n^  facilities 
corresponding  to  all  the  triplets  (i,  j,  k)  with  iel,  jeJ,  keK.  The 
transportation  costs  from  a  triplet  (i,  j,  k)  e  F  to  i,  j  and  k  are  equal  to 


0.  All  other  transportation  costs  are  equal  to  1.  The  capacities  are  equal 


1 


to  3,  the  demands  to  1  and  the  fixed  costs  to  e  <  -.  3DM  has  a  solution  if 


and  only  if  CPLP  has  a  solution  with  value  n  e. 


,B 


To  see  that  Zp  or  Zg  can  be  computed  in  polynomial  time,  note  that 


0  <  £.  d.  <  z.  s.  <  pn,  where,  here,  n  denotes  the  number  of  facilities, 
li  j  j 


Since  p  is  fixed,  a  dynamic  programming  algorithm  for  solving  the  knapsack 
problem  defined  by  (T)  is  polynomial.  The  result  now  follows  from  Grbschel, 
Lovasz  and  Schrijver  [15]. 

To  see  that  Z^  and  Zc  are  MP-hard,  note  that  the  transformation  of  3DM 
described  above  is  valid  even  if  we  relax  the  capacity  constraints  (C)  and  (T) 
completely.  3DM  has  a  solution  if  and  only  if  Z™^  has  the  value  ne.  This 
completes  the  proof  of  (3.4). 

Now  we  prove  (3.5).  From  an  instance  of  CPLP  where  the  demands  and 
capacities  are  taken  from  {1,2},  we  construct  a  graph  as  follows.  For  each 
customer  i  with  demand  d^  we  construct  d^  nodes,  for  each  facility  j  with 
capacity  Sj  we  construct  Sj  nodes  and  for  each  facility  j  with  capacity  1  we 
construct  an  additional  node.  The  corresponding  node  sets  are  denoted  by  LI,  V 
and  W  respectively.  The  graph  has  an  edge  with  cost  c^j  joining  each  node  of 
U  associated  with  customer  i  to  each  node  of  V  associated  with  facility  j. 
For  each  facility  j  such  that  Sj  =  2,  the  two  nodes  of  V  associated  with  j  are 


joined  by  an  edge  with  cost  -  f,. 


Finally,  each  node  of  V  associated  with  a 


■_  ■» 


\ 
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facility  with  Sj  =  1  is  joined  to  the  corresponding  node  of  W  by  an  edge  of 
cost  -  fj.  Now  CPLP  is  equivalent  to  solving  a  minimum  cost  matching  problem 
in  this  graph,  where  each  node  of  U  is  required  to  be  exactly  matched  and  each 
node  of  V  U  W  must  meet  at  most  one  edge  of  the  matching.  (The  closed 

facilities  correspond  to  the  nodes  of  V  that  are  matched  to  a  node  of  W  or  to 

another  node  of  V).  (3.5)  follows  from  the  fact  that  the  matching  problem  can 

be  solved  in  polynomial  time  (Edmonds[9] ) .  D 

It  is  interesting  to  note  that,  although  Z  can  be  computed  in  polynomial 

time  under  the  assumption  made  in  (3.5),  it  is  not  clear  whether  and  Z^ 

can.  In  fact,  there  are  values  of  the  Lagrange  multipliers  for  which  the 
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Lagrangian  relaxations  2^  and  2^  are  NP-hard  to  solve.  (For  example, 
the  reduction  of  3DM  used  for  proving  (3.4)  is  still  valid  when  v  =  0.) 

4.  Computational  Results 

I  t 

Our  computational  study  compares  the  relaxations  Z  ,  Z^,  Z^  and  Z^  to 
the  optimum  value  Z  of  CPLP.  We  also  compare  a  dual-based  heuristic  to  greedy 
and  interchange  heuristics.  One  objective  of  the  study  is  to  test  the 
tightness  of  the  bounds  as  a  function  of  parameters  of  the  problem.  In 
particular,  we  consider  the  influence  of  tight  capacity  constraints. 

Zq  is  computed  using  a  subgradient  algorithm  stopped  after  a  maximum  of 
200  iterations.  The  knapsack  subproblems  are  solved  using  the  algorithm  of 

Fayard  and  Plateau[11],  To  compute  Z*,  we  solve  the  Lagrangian  dual  giving 

T  IT 

ZDC  by  a  subgradient  algorithm  (recall  that  Z  =  ZD(,).  To  obtain 

T 

Z^,  Z  and  Z,  we  implement  the  cross  decomposition  of  Van  Roy  [30].  The 
optimum  value  Z  is  computed  using  a  branch-and-bound  algorithm  based  on  the 
bound  Z £.  The  bound  Zq  itself  is  calculated  by  solving  iteratively  a  sequence 


24 


of  problems  Z^u.  For  each  value  of  u,  Z^u  is  obtained  using  a  branch-and- 

T 

bound  algorithm  based  on  Z^,  .  Finally,  the  simple  plant  location  problem 
T 

that  yields  Z^  is  solved  using  the  dual  ascent  procedure  of  Erlenkotter 
[10]  which  also  involves  branch  and  bound.  Clearly,  due  to  the  triply  nested 
branch-and-bound  algorithms,  cross  decomposition  is  most  effective  when  the 
duality  gap  Z  -  Z£  is  small.  Dual  heuristics  are  needed  to  accelerate  the 

search  for  the  optimal  Lagrange  multipliers.  Lacking  this,  it  is  not  obvious 

how  to  generate  good  dual  information  from  the  primal  subproblems.  For  the 
harder  problem  instances,  the  cross  decomposition  method  often  reduces  to  a 
rather  inefficient  column  generation  method,  as  corroborated  by  Van  Roy. 

On  the  heuristic  side,  we  implemented  DROP-HI  and  VS M  -  a  greedy 
heuristic  starting  from  all  facilities  open  and  an  interchange  heuristic, 
respectively  -  as  described  in  Jacobsen[20j .  The  implementation  suggested  by 
Khumawala  [21]  gives  similar,  although  slightly  inferior,  results.  The 

computational  results  provided  below  follow  Jacobsen  [20].  We  also  consider  a 
dual-based  heuristic,  denoted  by  (H).  At  each  step  of  the  subgradient 
algorithm  for  computing  Zq,  a  solution  (x^,  y^)  of  the  Lagrangian  subprobiem 
is  found.  If  it  is  feasible  for  the  original  problem,  it  is  an  optimum 

solution.  Otherwise,  consider  y^  as  the  set  of  open  facilities  and  solve  a 
transportation  problem  to  find  a  feasible  solution  (x^,  y^).  The  best  such 
solution  found  during  the  subgradient  iterations  is  taken  as  the  solution  of 
heuristic  (H).  We  have  observed  that,  typically,  only  few  different  vectors 
y^  are  generated  during  the  algorithm  and  so  only  few  transportation  problems 
have  to  be  solved.  We  use  the  code  of  Srinivasan  and  Thompson  [29]  for 
solving  the  transportation  problems. 

To  validate  our  algorithms,  we  solve  the  classical  test  problems  of  Kuehn 
and  Hamburger  [23]  and  are  able  to  confirm  the  computational  results  already 


reported  in  the  literature.  For  example  our  evaluation  of  — —  coincides 

with  e  in  table  II  of  Van  Roy  [30]  in  all  but  one  of  the  instances.  (In 

problem  24,  we  found  .18$  instead  of  .21$  but  Van  Roy  truncated  his  branch- 

and-bound  after  2000  nodes  whereas  we  did  not.)  Van  Roy  does  not  report 

T  T 

Z  .  WeTfound  that  Z„  =  Z  for  all  the  problems  except  problems  11  and  12, 

L~Lz 

where  — —  =  .074$  and  .277$  respectively.  Regarding  Zp,  our  results  are 
better  than  expected.  The  reason  is  that,  although  Nauss  [25]  recommends  to 
compute  the  bound  Z^  in  his  paper,  he  reports  computational  results  that  are 
not  based  on  Zq.  More  specifically,  Nauss  relaxed  the  integrality 

requirements  in  the  knapsack  subproblems  needed  to  compute  Z^  (see  first 
sentence  of  second  paragraph  in  section  COMPUTATIONAL  RESULTS  of  [25]).  The 
reason  for  doing  this  is  unclear,  since  the  knapsack  subproblems  are  very  easy 
to  solve  when  all  the  Sj's  are  equal.  As  a  result,  Nauss  computed  Z1  instead 
of  Zq.  Indeed,  in  Table  I  of  [25],  the  gaps  closed  by  LGR1  and  LGR2  after 
subgradient  optimization  are  very  similar.  (By  Theorem  1,  they  should 
eventually  be  identical,  after  sufficiently  many  steps  of  the  subgradient 
algorithm.)  By  solving  true  knapsack  subproblems,  we  find  that  Zn  often 


greatly  improves  on 


For  example,  Zn  =  Z  for  problems  6,  8  and  9  as 


numbered  by  Nauss.  Going  back  to  the  numbering  of  Van  Roy,  we  find  that  Z^ 

=  Z  for  problems  1-6.  For  the  more  tightly  capacitated  problems  7-12,  we  find 
Z-ZD 

— =  .06$,  .20$,  .64$,  .95$,  .66$  and  .36$  respectively.  For  these  12 

instances,  the  heuristic  (H)  always  finds  the  optimum  solution.  The  heuristics 
DROP-HI  and  VSM  perform  reasonably  well  (see  Jacobsen  [20]). 

To  complement  these  results,  we  run  the  algorithms  on  a  different  and 
larger  set  of  test  problems. 

We  generate  5  sets  of  30  problems.  In  each  set,  there  are  5  problems  of 


each  of  the  following  sizes  -  25x8,  25x16,  25x25,  50x16,  50x33,  50x50  where 


the  first  number  is  the  number  of  clients  and  the  second  is  the  number  of 

facilities.  The  transportation  costs  are  computed  by  generating  points 

representing  the  clients  and  facilities  uniformly  in  a  unit  square.  The 

Euclidean  distance  between  them  is  multiplied  by  10  to  define  unit 

transportation  costs.  The  demands  are  generated  from  U[5,35],  where  U[a,b] 

denotes  the  uniform  distribution  in  the  interval  [a,b).  The  capacities  Sj  are 

generated  from  U [ 10, 1 60 ]  and  the  fixed  costs  by  the  formula  fj  =  U [ 0 , 90  ]  + 

U[  100, 110]  /Sj  to  reflect  economies  of  scale.  From  the  above  data,  5 

different  sets  of  problems  are  generated  by  scaling  the  capacities  using  the 
z  s . 

ratios  =  1.5,  2,  3,  5,  10.  Finally  f^  is  doubled  for  the  first  2 

l  ui 

sets.  Although  the  above  choices  are  somewhat  arbitrary,  our  objective  is  to 
generate  a  reasonably  difficult  set  of  test  problems.  The  actual  data  are 
available  from  the  authors. 

We  report  our  computational  experience  with  these  5  sets  of  problems  in 
Tables  I-V.  In  each  case  we  give  the  relative  error  of  the  bound  to  the 

z-zc 

optimum  value.  For  a  relaxation  such  as  Zq  we  report  — — ,  and  for  a 

Z ( h ) -Z 

heuristic  with  value  Z(H)  we  report  —  ; . 


Table  I  Ej  sj  =  1-5  £ 

.  d. 

1  1 

Relaxations 

Heuristics 

Relative  Error  % 

ZC 

7  ZT  7  ^ 

H 

VSM  DR0P-HI 

0 

0 

0 

0 

0 


0.37 

0 

0 

0.12 

0.01 


0 

0 

0.04 

0 

0.07 


0 

0 

0 

0.26 

0 


0.  10 
0.05 
0 

0.  12 


0 

5.73 

6.22 

13.13 

8.82 


15.69 

6.23 

13.30 

10.92 

0.29 


5.95 

10.02 

9.18 

15.01 

12.57 


5.30 

7.75 

4.08 

2.77 

0.48 


25.52 

19.35 

18.58 

5.77 


17.27 
5.73 
11.14 
13.  13 
8.82 


13.30 

15.20 

5.36 


30.01 

20.96 

26.96 

38.25 

26.36 


17.97 

23.71 

24.42 

16.78 

22.86 

26.24 

24.10 

28.02 

26.10 


0.05 

1 .61 

2.89 

0 

11.72 

29.98 

0 

0.06 

0.25 

1 .45 

0 

# 

27.63 

0 

0.07 

0.09 

19.62 

0 

* 

26.08 

0 

0.02 

0.38 

32.96 

0 

* 

24.44 

0.01 

0.02 

0.19 

0.41 

0 

# 

28.88 

0 

0. 10 

0 

11.80 

0 

# 

27.58 

0.04 

0.07 

1.97 

9.83 

0 

10.57 

21.46 

50  x  50 


Average  % 


time  limit  of  120  seconds  on  DEC  20  exceeded. 
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Table  II  I .  s ,  =  2  E.  d. 

J  J  11 


Relaxations 

Heuristics 

Relative  Error  % 

ZC 

ZD 

ZT 

ZC 

Z1 

H 

VSM 

DROP-HI 

0 

0 

6.22 

0 

0 

24.05 

0.75 

0.79 

■  -  • 

10.32 

0 

0.82 

8.40 

25  x  8 

0 

0.05 

7.80 

0 

4.27 

14.40 

0 

0.02 

2.25 

0 

0 

35.71 

0 

0 

5.26 

0 

19.54 

21.84 

0.19 

4.91 

4.97 

0 

10.36 

41.98 

0.09 

1.43 

1 .47 

0 

22.36 

25  x  16 

0 

0 

2.09 

2.27 

0 

8.75 

38.64 

0.15 

0. 15 

4.36 

4.36 

0 

14.30 

20.91 

0 

0 

2.74 

3.34 

0 

18.59 

44.04 

0.05 

0.12 

2.28 

2.28 

0 

4.82 

34.68 

0.18 

0.23 

1 .42 

1.57 

0 

9.56 

29.89 

25  x  25 

0 

0.04 

0.78 

0.84 

0 

20.25 

40.33 

0 

0.03 

1.76 

2.38 

0 

13.30 

37.25 

0.1 1 

0. 12 

3.04 

3.07 

0 

6.76 

41.36 

0 

0. 19 

3 . 66 

3.89 

0 

13.  19 

27.92 

0.01 

0.17 

2.55 

2.57 

0 

5.93 

23.1 1 

50  x  16 

0 

0 

0.46 

0.49 

0 

11.05 

25.08 

0.42 

0.43 

2.45 

11.85 

0 

14.  12 

14.12 

0 

0 

0 

0 

0 

7.61 

41.27 

0.04 

0.04 

1.50 

1 .50 

0 

10.84 

46 . 67 

0.09 

0.09 

1.15 

1.19 

0 

1 1  .81 

32.34 

50  x  33 

0.08 

0.  10 

1.52 

1.86 

0 

19.98 

47.39 

0.01 

0.02 

0.67 

0.81 

0 

14.38 

28.48 

0.09 

0.09 

1.34 

1 .34 

0 

12.82 

39.08 

0 

0 

0.23 

0.34 

0 

# 

36.27 

0.21 

0.27 

1 .01 

1 .25 

0.01 

# 

41.60 

50  x  50 

0 

0.02 

0.14 

1.43 

0 

* 

45.38 

0 

0.03 

0.13 

0.36 

0 

# 

39.02 

0 

0.04 

0.02 

0. 13 

0 

# 

42.78 

Average  % 

0.08 

0.1 1 

2.43 

2.91 

0 

10.64 

33.21 

29 


30 


Table  IV  z  s  =  5  Z.  d. 


Relaxations 
!  „ 


Heuristics 


zc 

7 

ZT 

Lz 

■■ 

H 

VSM 

DROP-HI 

0 

0 

11.61 

11.75 

0 

0 

3.92 

0 

0 

5.95 

5.95 

0 

3.98 

10.64 

0 

0.31 

5.39 

5.65 

0 

0 

18.88 

0 

2.28 

11.40 

11.66 

0 

0 

11.90 

0 

0 

0 

0 

0 

0 

23.92 

0 

0 

1.57 

2.07 

0 

0 

59.50 

0 

0 

0.22 

0.82 

0 

1.61 

49.81 

0 

0 

3.72 

4.06 

0 

4.55 

37.36 

0 

0 

1.76 

2.01 

0 

3.18 

25.06 

0.12 

0.12 

1.54 

2.34 

0 

0 

40.53 

0 

1 .07 

3.76 

3.89 

0 

3.33 

52.94 

0.19 

0.74 

2.68 

12.10 

0 

13.31 

32.80 

0.27 

2.57 

5.90 

5.90 

0.07 

7.11 

59.22 

0.30 

0.53 

2.43 

2.90 

0 

0 

64.85 

0 

0.22 

2.11 

2.20 

0 

2.22 

34.07 

0 

0 

0.51 

1  .40 

0 

0 

30.36 

0 

0.70 

3.27 

3.68 

0 

2.37 

28.91 

0.38 

0.38 

1 .87 

3.92 

0 

0 

31.11 

0.65 

2.10 

4.73 

4.90 

0 

0 

22.01 

0 

0.85 

4.13 

4.31 

0 

1.76 

25.55 

0.59  3.04 


42.37 


Average  % 


77. 

.06 

81, 

.07 

76. 

.  16 

75. 

.74 

63. 

.  1 1 

42. 

.81 

Some  remarks  can  be  made  about  Tables  I-V.  The  relaxations  Zq  and  Zp  are 
both  significantly  stronger  than  and  Z  for  this  set  of  test  problems.  Zq 
is  about  twice  as  strong  as  Zq  for  the  tightly  capacitated  problems  and  this 
advantage  increases  for  the  less  capacitated  ones.  The  quality  of  both  bounds 


deteriorates  in  the  intermediate  range,  i.e.  I^s^/I^d^  =  3  and  5,  but  Z^ 


S  V  *v  *.  •/  ^  »>w  V 
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never  exceeds  \%  and  ZD  never  exceeds  3>%.  Whereas  for  the  Kuehn  and  Hamburger 
test  problems  Z£  is  frequently  equal  to  Z^,  here  Z£  is  typically  much 
weaker.  For  the  tightly  capacitated  problems,  the  negative  effect  of  removing 
the  knapsack  constraint  (T)  is  more  pronounced  when  there  are  few  facilities 
(e.g.  8  or  16).  Our  subgradient  algorithm  for  computing  had  difficulties 
converging  for  the  tightly  capacitated  problems  of  Table  I.  So  the  large  gaps 
appearing  there  should  be  interpreted  with  caution.  DROP-HI  and  VSM  are  the 
two  heuristics  recommended  by  Jacobsen,  based  on  tests  with  the  Kuehn  and 
Hamburger  problems.  Their  performance  is  very  different  when  the  facilities 
do  not  have  identical  capacities  and  fixed  costs,  as  pointed  out  by  Domschke 
and  Drexl  [8].  Our  own  computational  experience  serves  to  reiterate  this 
warning.  DROP-HI  tends  to  close  the  large  facilities  first,  clearly  not  a 
good  strategy.  VSM  performed  reasonably  well  on  the  less  capacitated 
problems.  However,  the  clear  winner  is  heuristic  (H),  failing  to  find  the 
optimum  solution  in  only  10  out  of  150  cases.  This  is  particularly  remarkable 
considering  that  the  time  needed  to  compute  Z^  and  to  find  the  heuristic 
solution  (H)  only  grows  moderately  with  problem  size.  By  contrast,  VSM 
exceeded  the  time  limit  of  120  seconds  on  a  DEC  20  in  thirteen  instances. 
Computing  times  are  provided  in  Table  VI.  More  extensive  computational 
results  are  reported  in  Sridharan  [28], 


V 


Table  VI  Computing  Times  (Average  of  25  problems). 


Problem  Size 


Seconds  on  IBM  3081 
,  „T 


Seconds  on  DEC  20 


ZD  and  ZH 


DROP-HI 


25  X  8 

0.9 

1.0 

2.  1 

,.7 

0.6 

0.  1 

25  X  16 

14 

12 

5.7 

3.9 

2.8 

0.2 

25  X  25 

140 

69 

1 1 

9.0 

9.6 

0.5 

50  X  16 

34 

23 

8.7 

9.6 

5.7 

0.5 

50  X  33 

1900 

570 

28 

23 

41 

1.7 

50  X  50 

6800 

5200 

62 

37 

>92* 

3.7 

3  instances  exceeded  the  time  limit  of  120  seco"dc  on  a  DEC  20. 


Probabilistic  Results 
In  this  section  we  assume  that: 


)  n  clients  are  distributed  independently  and  uniformly  in  the  unit 
square, 


i)  every  client  location  is  also  a  potential  site  for  locating  a  facility, 


ii)  Cy  is  proportional  to  the  Euclidean  distance  between  points  i  and  j, 


v)  all  the  clients  have  the  same  demand,  independent  of  n, 


)  all  the  facilities  have  the  same  fixed  cost  f  =  f(n)  and  the  same 


capacity  s  =  s(n) . 


i 


{jr^-y'^vrr.Typr.-v^- 
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(vi)  s/log  n  -*■<=  as  n-*<= , 


(vii)  ne-^^  <f<  rJ-e  for  some  fixed  e>0. 


Although  these  assumptions  are  restrictive,  we  believe  that  they 
represent  an  interesting  starting  point  for  a  probabilistic  comparison  of 
heuristics  and  relaxations  for  the  CPLP.  In  fact,  the  results  we  are  about  to 
present  still  hold  under  weaker  assumptions.  In  (i),  the  shape  of  the  region 
can  be  almost  anything  (the  neighborhood  of  the  boundary  must  be  of  measure  0) 
and  the  distribution  can  be  nonuniform  as  long  as  mass  is  not  concentrated. 
In  (ii),  the  potential  sites  need  not  coincide  with  the  clients  as  long  as 
they  have  a  similar  distribution.  In  (iii),  the  Euclidean  distance  can  be 
computed  after  a  linear  transformation  of  the  region.  In  (iv),  the  demands 
can  be  different  for  each  client  as  long  as  they  are  bounded  independently  of 
n.  In  (v),  the  fixed  values  f  and  s  can  each  be  replaced  by  a  small  range  of 
values.  The  assumption  in  (vii)  can  be  slightly  relaxed  to  u>//n  <f<  n/w 
where  as  n-“.  Our  goal  in  this  section  is  only  to  give  a  flavor  of 

typical  results  and  assumptions  (i)-(vii)  serve  this  purpose. 

Without  loss  of  generality,  we  take  the  constant  of  proportionality  in 
(iii)  to  be  1  and  the  demand  in  (iv)  to  be  1.  Then,  given  f  and  s  considered 
as  functions  of  n,  the  value  Z  =  Z(n)  is  a  random  variable  for  every  n  >  1. 
Similarly,  relaxations  or  heuristics  yield  random  variables  Z^(n)  or  Z^(n)  for 
every  n  >  1.  In  this  section  we  study  the  probabilistic  convergence  of  such 
sequences  of  random  variables  as  n  -> 

First,  we  introduce  some  notation.  Let  W(n)  be  a  sequence  of  random 
variables  for  n  >  1.  We  say  that  W(n)  ->  w  almost  surely  if  and  only  if 


(5.1)  Pr[lim  sup  W(n)  =  w  =  lim  inf  W(n)]  =  1. 

If  we  only  have  Pr[lim  sup  W(n)  <w]  =  1,  we  say  that  W(n)  <  w  almost 
surely.  A  sufficient  condition  for  (5.1)  is  ( Borel-Cantelli  lemma) 

(5.2)  1^]  Pr[|W(n)-w|  >  e]  <  ”,  for  every  fixed  e>0. 

Let  V(n)  and  W(n)  be  two  sequences  of  random  variables  for  n  >  1.  We 
write  W(n)  -  V(n)  almost  surely  if  and  only  if  W(n)/V(n)  ->  1  almost 
surely.  Frequently,  we  simply  write  W  ~  V  almost  surely. 

In  this  section,  we  show  thac  the  relaxations  of  the  CPLP  introduced 
earlier  fall  into  3  classes  according  to  the  probabilistic  convergence  of 
(Z-ZR)/Z.  Let  a  =  .32887076 _ 

Theorem  5 

Assume  (i)  -  (vii).  Then 

T 

(a)  Z  -  Z  almost  surely, 

(b)  ~  Z*  almost  surely  and 
ZJ/Z  >  .99810744  almost  surely, 

(c)  Z®1  =  fn/s  and 


Z®*/Z  <  .33333334  almost  surely  for  f  <  (asl^^/Zn. 


For  f  >  (as)  //n,  the  ratio  Z  /Z  improves  above  33%  but  the  bound 

or 

Z  is  still  significantly  weaker  than  the  ones  considered  in  (a)  and  (b). 
The  3  parts  of  the  theorem  are  proved  in  Lemmas  6,  7  and  8  respectively.  Our 
analysis  uses  the  following  results. 

Lemma  3:  (Hoeffding  [19])  If  Y^,  ...,  Yn  are  independent  random  variables  and 
0  <  <  1  for  i=1,...,n,  then,  for  0  <  e  <  1, 

Prt^Y.  >  ( 1+e )n ]  <  e-*2">*/3 

where  p  is  the  expected  value  of  (^n^Y^)/n. 

Papadimitriou  [26]  introduced  the  following  K-honeycomb  heuristic. 
Partition  the  plane  into  a  regular  hexagonal  tiling,  where  each  hexagon 
has  the  same  area  1/K.  For  each  hexagon  that  intersects  the  set  of 
clients,  let  P^  be  the  client  closest  to  its  center.  Choose  this  set  [P^]  as 
the  set  of  open  facilities.  We  denote  by  Z-^  the  value  of  this  heuristic. 

Papadimitriou  [26]  proved  that  Z^  is  near-optimai  for  the  K-median 
problem.  A  similar  result  holds  for  the  SPLP  for  an  appropriate  choice  of  K, 
as  shown  in  Ahn,  Cooper,  Cornuejols  and  Frieze  [1], 

Lemma  4 :  Assume  (i)  -  (vii)  and  s  =  =°.  Then 

(a)  Z^  ~  Z  almost  surely  for  K  ~  a(n/f)273, 

(b)  Z  -  3a(fn2)173 


almost  surely. 


To  approximate  the  linear  programming  relaxation  of  SPLP,  we  define  the 
following  feasible  fractional  solution  for  each  K=1,...,n.  For  j=1,...,n,  set 
=  K/n.  For  i=1,...,n,  define  j*  to  be  a  facility  such  that  c^  <  c^*  for 
at  most  n/K  values  of  j  and  c„  <  c.^*  for  at  least  n/K  +  1  values  of  j; 
take  xij  =  K/n  for  n/K  facilities  j  *  j*  such  that  c^  <  c^*,  take  x^#  =  1 
n/K  K/n  and  x  —  =  0  for  every  other  value  of  j .  We  call  this  solution  the 
K/n-uniform  fractional  solution  and  denote  by  Zp  =  2^1.  c^x_  +  Kf  its 
value . 

Lemma  5:  (Ahn,  Cooper,  Cornuejols  and  Frieze  [1])  Assume  (i)  -  (vii)  and 
s=”.  Then 

(a)  Zc  -  Z1  almost  surely  for  K  -  s(n/f)2/j, 

F 

(b)  Z1  ~  38(fn2)1/'^  almost  surely, 

where  8  =  .32824836  ...  =  1/(9^) 

T 

Lemma  6:  Z  -  Z  almost  surely. 

T 

Proof:  First,  we  find  a  lower  bound  on  Zq. 

zl  =  Max  Min  Z . Z  . ( c .  ,+v  )x .  +  Z  .(f-sv  )y 

c  0  i  J  ij  J  ij  J  J  J 

J  Z . x .  =  1  for  all  j 

i  ij 

0  <  x .  .  <  y  .  for  al  i ,  j 
ij  J 

y,  e  (0,1} 


for  all  j. 


Setting  Vj  =  v  for  all  j,  we  get,  for  every  v  >  0, 


(5.3)  zl  >  nv  +  Z(v)  where 


Min  z.z.c..x..  +  r.(f-sv)y. 
i  J  ij  ij  J  J 


T.  .X  .  .=  1 

ij 

0  <  x .  .  <  y  . 
ij  J 

yj  e  {0,1} 


for  all  j 
for  all  i,  j 


for  all  j. 


This  is  an  SPLP.  Let  e  ->  0  as  n  ->  =>.  The  actual  value  of  e  will  be  given 
later.  By  Lemma  4(a),  the  K-honeycomb  heuristic  satisfies 


(5.4)  Zv  ~  Z(v)  almost  surely  for  K  and  v  defineo  as  follows 


K  =  n(1+E)/s  =  a(n/( f-sv) )^/^  when  a(n/f)^/^  <  n(1+e)/s 


K  =  a(n/f)^/3  an(j  v  =  0  when  a(n/f)^'^  >  n(1  +  E)/s. 


We  claim  that  the  K-honeycomb  heuristic  is  almost  surely  feasible  to  CPLP 
for  this  choice  of  K  and  v.  To  prove  this,  it  suffices  to  show  that  there  are 
almost  surely  at  most  s  clients  in  each  hexagon.  The  number  N  of  clients  in 
an  hexagon  of  area  1/K  has  the  binomial  distribution  8(n;  1/K).  By  Lemma  3, 


Pr[N  >  (  1  +  e )n/K ]  <  e 


-e  n/3K 


Let  E  be  the  event  that  at  least  one  hexagon  of  the  tiling  contains  more  than 
s  clients.  Since  K  ->  ”,  the  number  of  hexagons  of  the  tiling  that  intersect 
the  unit  square  is  (1+o(1))K.  Thus 

2 

Pr [E j  <  ( 1 +o( 1 ) )Ke_e  n/3K. 

p 

Let  us  choose  e  so  that  e  n/3K  =  31ogn.  It  is  easy  to  check  that  e-<-0  as  n-*-“. 
So 

Pr[E]  <  (1  +  o(1))K/n^  <  1/n^. 

This  shows  that  the  K-honeycomb  heuristic  is  almost  surely  feasible  for  CPLP 
for  the  above  choice  of  K  and  v.  This  implies 

(5.5)  Z  <  +  Ksv  almost  surely. 

The  inequalities  (5.3),  (5.4),  (5.5)  imply 

Zq  >  nv  +  Z^,  (1-o(1))  >  Z ( 1  — o (  1 )  +  nv  -  Ksv(1-o(1)  almost  surely. 

Since  n  =  Ks(1-o(1))  and  nv  <  Z,  we  get 

T 

Zq  >  Z  ( 1  — o (  1 )  )  aimost  surely. 

Since  zj  is  a  relaxation  of  CPLP,  z£  <  Z.  So 


T 

Z^  ~  Z  almost  surely. 


□ 


Lemma  7 : 


l} /Z  ->  .99810744  ...  almost  surely  when  f  <  s^^/3/iTn 
Z^/Z  >  .99810744  almost  surely  when  f  >  s^^/3/nn 
Z^  -  2}  almost  surely. 


Proof : 


min  lie .  .x .  .  +  Ify  . 

iJ  iJ  J 

r .x.  .  =  1 
J  ij 

I . x .  .  <  sy  . 

1  ij  J 

0  <  x .  .  <  y  .  <  1 
ij  J 


Consider  the  K/n-uniform  fractional  solution  of  Lemma  5,  for 
K  =  6(  n/f  )~>^  (  1  +  e  )  where  e  ->  0. 

Assuming  f  <  s-^^^/3/n,  we  get  n/K  <  s(1  +  e). 

So  this  solution  satisfies  the  capacity  constraints  and  therefore  is 
feasible  to  l} ,  i.e.  Z^  <  Zp. 

By  Lemma  5,  Z„  -  Z ^  ~  36(fn^)^^  almost  surely. 

Since  ZC*  <  Z1,  we  get 
Z1  -  BBCfn^)1^  almost  surely. 

The  inequality  f  <  s^/^/3/^n  implies  a(n/f)^^^  >  n(1+e)/s.  We  have  shown 
in  Lemma  6  that,  under  this  condition,  the  K-honeycomb  heuristic  is  almost 
surely  feasible  to  CPLP  for  K  =  a(n/f)^/3 


By  Lemma4,  this  implies 


Z  -  3a(fn  ) 


almost  surely. 


Thus,  when  f  <  s^^^/3/Ttn,  7} /Z  ~  &/a  almost  surely. 

Since  8/a  =  .99810744  ...,  we  have  shown  the  first  part  of  Lemma  7. 
When  f  >  s^/^/3/^n,  note  that  we  can  write 


„  min  lie.  .x.  .  +  I(f-sv)y. 
/>0  ij  ij 

I  ,x .  .  -  1 
J  ij 

I . x .  .  <  sy  . 

1  ij  J 

0  <  x.  .  <  y  .  <  1 . 
ij  J 


Now  consider  the  K/n-uniform  fractional  solution  for  K 
(n/(  f-sv)  )^/'^(  1+e )  where  e  ->  0.  By  Lemma  5,  its  value  satisfies: 


Zp  -  3 B [ (f-sv)n2l ^  almost  surely. 


Therefore 


I  2  1/3 

Z  >  maxv>Q3B[ ( f-sv)n  ]  +■  nv  almost  surely. 


This  function  of  v  has  a  unique  maximum,  obtained  when: 


(5.6)  n/s  =  6(n/(f-sv) )^/^. 

This  implies  n/K  =  s/(1+e)  <  s.  So  the  above  K/n-uniform  fractional  solution 
is  feasible  for  Z*.  This  shows  that 


Z  ~  38[(f-sv)n  ]  +  nv  almost  surely 


where  v  is  defined  by  equation  (5.6). 


Similarly 


max  min  ZZc.  .x.  .  +  Z(f-sv)y  .  +  nv 
v>0  ij  ij 

Z .x.  .  =  1 
J  ij 

Z . x .  .  <  sy  . 
i  ij  J 

0  <  x .  .  <  y  . 

J-J  J 

yj  £  (0,1). 


By  Lemma  4,  the  K-honeycomb  heuristic,  for  K  =  a(n/( f-sv) )2/ 1+e )  and  t 


satisfies 


ZK  ~  3a[ ( f-sv )n^ ] 1 ^  almost  surely. 


Z  >  max  >n3a[ ( f-sv)n2 ] 1 +  nv  almost  surely. 


The  value  of  v  that  yields  the  maximum  is  given  by 


n/s  =  a(n/(  f-sv)  )2//^ 


.•  V  /  -v  >  .*  -*  V  -*  / 


This  implies  n/K  =  s/(1+e).  By  an  appropriate  choice  of  e  we  can  show  as 
in  Lemma  6  that  the  K-honeycomb  heuristic  is  almost  surely  feasible.  Thus 


2,1/3 


Z  ~  maxv>Q3a[  ( f-sv)n  ]  +■  nv  almost  surely. 


Therefore 


Z  <  3o[ ( f-sv)n2 ] 1 ^  +  nv  almost  surely 


where  v  is  the  value  defined  by  equation  (5.6) 


So  Z*VZ  >  (38[  (f-sv)n^]  ^/^+nv)/(3a[(f-sv)n^]  ^'3+nv)  >  8/a  almost  surely 


since  8  <  a. 


Finally  we  prove  Z~  ~  Z  .  By  (2.3),  Zp  =  ZpC> 


Z“  =  min  ZIc . .x. .  +  Ify . 
DC  ij  ij 


( x ,y )  e  P(DC)nP( I NT) . 


Zqq  reduces  to  a  linear  program  very  similar  to  Z.  The  only  difference 
is  the  constraint 


(5.7)  Z  y j  >  n/s  in  Z^q  instead  of  l  y^  >  n/s  in  Z*. 


Consider  a  K/n-uniform  fractional  solution. 

When  e(n/f)^3  >  n(1+o(1))/s,  set  K  ~  B(n/f)^/^.  Then  neither  of  the 
constraints  (5.7)  is  binding.  So,  by  Lemma  5, 


Zp  ~  Z  ~  Z^  almost  surely. 


When  8(n/f)2//3  <  n(1+o(1))/s,  then  set  K  =  [n/s].  Note  that  K  ~  n/s  as 
n/s  ->  <*>. 

So  Zp  ~  -  Zp  almost  surely  by  the  analysis  made  earlier  to  compute 

Z*  when  f  >  s^^/S/im.  Q 


Lemma  8 


ZBI  =  fn/s 


ZBI/Z  <  .3333334  almost  surely  for  f  <  (as)^/2//n. 


Proof :  The  constraint  lyj  >  n/s  is  feasible  for  Z  .  So 


ZBi  =  ZZ  c,  : x j  ^  +  fry  .  >  fn/s. 


Now  consider  the  following  solution. 


yj  =  1/s  for  all  j,  x^  =  1,  x^j  =  0  for  i  *  j.  It  is  feasible  for  Zf 


and  its  value  is  f  n/s.  So  ZB^  =  f  n/s. 


When  f  <  (as)^/'^//n ,  then  Z  -  3a(fn‘^)'''/3  almost  surely. 


Z  >  3(1  —0  (  1 ) )  (f2n)  ^^(fn2)  ^3/s  =  3(  1 -o(  1 ) )  fn/s  almost  surely. 


This  shows  ZB*/Z  <  1/3+o(1)  almost  surely. 


6.  Conclusion 


In  this  paper  we  have  studied  relaxations  of  CPLP  from  three  different 
angles:  (i)  inequalities  among  the  corresponding  bounds,  (ii)  computational 
experiments  and  (iii)  probabilistic  analysis.  The  results  in  (i)  can  be 
thought  of  as  a  ready  reckoner  to  indicate  the  relaxation  that  is  easiest  to 
solve  among  those  that  yield  the  same  bound.  We  also  showed  that  relaxations 
based  on  the  variable  splitting  method  do  not  provide  different  bounds  than 
the  classical  Lagrangian  relaxations. 

In  the  probabilistic  study,  where  we  assume  equal  fixed  costs  and 
capacities  at  all  facilities,  we  showed  that  the  relaxations  fall  into  three 
categories  according  to  their  strength.  The  weak  LP  relaxation  can  be  solved 
analytically  but  the  bound  is  indeed  very  weak.  The  strong  LP  relaxation 
almost  surely  gives  a  bound  within  one  fifth  of  one  percent  of  the  optimum. 
The  relaxation  proposed  by  Nauss  falls  into  the  same  category.  These 
relations  are  relatively  easy  to  solve.  Finally,  the  relaxation  proposed  by 
Van  Roy  is  extremely  tight  as  it  almost  surely  provides  a  bound  that  goes  to  Z 
as  the  size  of  the  problem  increases.  But  this  relaxation  is  strongly  NP-hard 
to  compute. 

The  computational  results  complement  our  theoretical  results.  For  and 
Zq  we  found  that  problems  that  are  either  loosely  or  tightly  capacitated  have 
smaller  duality  gaps  compared  to  those  in  the  intermediate  range.  Based  on 
our  computational  results  we  suggest  the  use  of  a  relaxation  based  heuristic 
(H)  to  solve  CPLP.  It  provides  considerably  better  solutions  that  the  VSM 
heuristic  and  is  also  superior  on  computational  time,  at  least  on  the  large 
problems.  Further,  it  is  obtained  in  conjunction  with  a  lower  bound. 
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